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Abstract 

We study nematic liquid crystal configurations in a prototype bistable 
device - the Post Aligned Bistable Nematic (PABN) cell [9]. Work- 
ing within the Oseen-Frank continuum model, we describe the liquid 
crystal configuration by a unit-vector field n, in a model version of the 
PABN cell. Firstly, we identify four distinct topologies in this geome- 
try. We explicitly construct trial configurations with these topologies 
which are used as initial conditions for a numerical solver, based on the 
finite-element method. The morphologies and energetics of the cor- 
responding numerical solutions qualitatively agree with experimental 
observations and suggest a topological mechanism for bistability in 
the PABN cell geometry. 
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1 Introduction. 



Liquid crystals are an intermediate phase of matter between the solid and 
liquid states. In the simplest liquid crystal phase, the nematic phase, the 
constituent rod-like molecules tend to align along a locally preferred direc- 
tion. This mean direction of molecular alignment is described by a director 
field n(r), which is an unoriented unit- vector field so that the sign of n has 
no physical significance pp. 

The existence of a locally preferred direction and the resulting anisotropic 
optical properties make liquid crystals very suitable for display devices. Most 
of the liquid crystal displays (LCDs) in use today, such as the Twisted Ne- 
matic and the Super Twisted Nematic, are monostable [3J. They can sup- 
port two optically contrasting states, only one of which is stable without an 
applied field. Recently, there has been considerable interest in developing 
bistable display technologies, where there are two or more stable, optically 
contrasting states [!],[§]■ Here, power is needed only to switch between the 
different states but not to maintain them. 

Bistable nematic LCDs typically use a combination of complex surface 
morphologies and surface treatments to stabilize the different states jU El [7] . 
This paper focuses on one such bistable device - the Post Aligned Bistable 
Nematic (PABN) device [§]• The PABN device consists of a liquid crystal 
layer sandwiched between two substrates. The lower substrate is featured 
by an array of microscopic posts, as shown in Figure (00). The boundary 
conditions are a mixture of tangent and normal conditions in various parts 
of the geometry (referred to as homogeneous and homeotropic respectively 
in liquid crystal literature pQ). Tangent boundary conditions on a surface 
constrain the director n to be in the plane of the surface whereas normal 
boundary conditions constrain n to be perpendicular to the surface. 

For a range of post heights, the PABN device is experimentally observed 
to be bistable [9j. It supports two optically contrasting states with long- 
term stability - one bright and the other dark when viewed between crossed 
polarizers. Optical modelling suggests that the dark state corresponds to a 
liquid crystal configuration that tilts strongly around the microscopic posts 
whereas the bright state corresponds to a suppressed tilt profile around the 
posts ie. a comparatively planar configuration [9]. The high-tilt configuration 
is referred to as the tilted state and the low-tilt configuration as the planar 
state in the remainder of the paper. 

We use topological arguments to study static director configurations in 
the PABN geometry. The admissible configurations in this geometry can be 
partitioned into distinct topological classes (for details see [H] and [13]) and 
one might expect a stable configuration (or local energy minimizer) for every 
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such class. The key point is that topologically distinct stable configurations 
cannot be continuously deformed into each other and switching proceeds 
via creation or annihilation of topological defects and boundary condition 
violation P,[7]. Therefore, such topologically distinct stable configurations, 
if they exist and have similar free energies, lead to bistability in prototype 
devices. 

Our study is related to basic questions about the existence and properties 
of multiple equilibrium configurations in complex geometries. For polyhedral 
geometries (such as a rectangular prism), we have addressed several general 
questions about the topological classification, energy estimates and regularity 
of equilibrium configurations in previous work [12], [14]. In this paper, we 
apply similar methods to a complex non-convex polyhedral geometry - the 
PABN geometry. 

The paper is organized as follows. In Section (2), we describe our sim- 
plified model of the PABN device, based on the Oseen-Frank continuum 
model. In Section (3), we identify four distinct (but simple) topologies in 
this model geometry. We construct trial configurations with these topologies 
in Section (4). These trial configurations are used as initial conditions for a 
numerical solver of the equilibrium Euler-Lagrange equations (T41). We obtain 
numerical solutions for each topological class. The energies and morphologies 
of these solutions are then correlated to the observed physical phenomena. 
The numerical results and experimental observations are in qualitative agree- 
ment. 

2 Model. 

We work with a simplified model version of the PABN geometry, neglecting 
factors such as post tilt and replacing the rounded edges and corners in 
Figure (CQ) by sharp features. Further, we only consider configurations which 
are periodic in the array so that it suffices to look at what happens around a 
single post. (This is a realistic assumption since a typical pixel in the PABN 
device consists of thousands of these posts.) 

Our model geometry is displayed in Figure (jl]) - a single rectangular post 
of fixed square cross-section (L p x L p ) and variable height, h, inside a cell of 
fixed dimensions, L c x L c x H. The post and cell cross-sectional parameters, 
L p and L c respectively, are chosen so that L c = 2L p whereas the cell height 
H is fixed to be H = 3L C . These choices roughly agree with the actual device 
parameters. 

We model the liquid crystal configuration, outside the rectangular post, 
by a unit- vector field n(r). We impose periodic conditions on the cell bound- 
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aries so that 

n (x, y, z) = n (x + L c , y, z) = n (x, y + L c , z) etc. (1) 

The constraint ([T]) allows us to extend our results to a periodic array of posts. 
Further, n is taken to satisfy tangent boundary conditions on the bottom 
substrate and post surfaces and normal conditions on the top substrate, 
consistent with the boundary conditions in the device. 

Working within the Oseen-Frank continuum model, the liquid crystal 
energy is given by 

E [n] = I w (n, Vn) dV (2) 
Jv 

where 

w (n, Vn) = K x (V • n) 2 + K 2 (n • V x n) 2 + K 3 (n x V x n) 

+ (K 2 + K 4 ) (tr(Vn) 2 - (V ■ n) 2 ) (3) 

and the Kj are material-dependent elastic constants pQ. The stable config- 
urations then correspond to local minimizers of ([2D, subject to the imposed 
boundary conditions. 

3 Topology. 

Our aim in this paper is to identify and analyse topologically distinct stable 
configurations in the PABN cell (ie. configurations that cannot be contin- 
uously deformed into each other). In this section, we first study the two- 
dimensional Zenithally Bistable Nematic (ZBN) cell as an illustrative ex- 
ample [H [5]. The ZBN cell has a comparatively simple two-dimensional 
geometry where the topology is characterized by a single quantity - the pla- 
nar winding number [8j. Here, we look at the experimentally observed states 
and identify their distinct topologies. Then, we study the three-dimensional 
PABN cell, where the geometry is far more complex and the topology richer. 
We identify four different topological classes in this geometry. These classes 
generate tilted and comparatively planar profiles around the rectangular post 
and, hence, serve as good candidates for the topologies of the experimentally 
observed tilted and planar states. 

3.1 ZBN cell. 

The ZBN cell, like most liquid crystal cells, consists of a liquid crystal layer 
sandwiched between two substrates. The bottom substrate is planar whereas 
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the upper substrate is a mono-grating. Both substrates are treated to be 
homeotropic so that n is constrained to be normal to these surfaces. 

The ZBN cell supports two stable, optically contrasting configurations - 
the high-tilt almost vertical state which is dark when viewed under crossed 
polarizers and the low-tilt state, supporting greater bulk distortion, that is 
bright under crossed polarizers [HJ H] (tilt is measured with respect to the 
horizontal direction). The cell cross-section (with a wedge-shaped upper 
substrate) and the director profiles for the two observed states are shown in 
Figures (Eh) and (Gb). 

The topology of these states is characterized by the rotation of n along 
a path connecting the two sides of the grating, as shown in Figure ([2]). The 
high-tilt state has the minimum allowed rotation consistent with the bound- 
ary conditions whereas the low-tilt state exhibits greater intermediate rota- 
tion. In fact, the net rotations for the two states differ by tt or a half- winding 
number [8] . This then implies that the two states have different winding num- 
bers and are thus topologically distinct, leading to their long-term stabilities. 

It should be noted that n can, in principle, rotate by arbitrarily large 
amounts around the relief grating (between the fixed orientations at the sides) 
but such states have much higher energies and are unlikely to be observed in 
practice. 

3.2 PABN cell. 

Next, we consider the three-dimensional model PABN cell, in Figure (HI). 
Here, the topology cannot be merely characterized by the planar winding 
number, as in the ZBN case. In fact, there are at least three separate topo- 
logical invariants on the post edges, post faces and around the post vertices. 
In this section, we do not give a detailed description of the topological classifi- 
cation but identify four distinct, low-energy topologies. A systematic account 
is given in [TT] and [13] . 

We first recall that the tangent boundary conditions on the bottom sub- 
strate and the post surfaces imply that on these surfaces, n takes values 
tangent to these surfaces. Therefore, on the post edges, n is parallel to the 
edges and can take one of two possible values. The value of n, on a post edge, 
is defined to be the corresponding edge orientation [TT] . n is necessarily dis- 
continuous at the vertices, where three or more edges meet. For simplicity, 
we only consider configurations which are continuous everywhere away from 
the sharp post vertices. 

We first look at the top face of the post and the corresponding four 
horizontal edges. Up to symmetry and the sign of n, there are three distinct 
choices of the horizontal edge orientations on this face, as shown in Figure ([5]). 
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Figure 1: The PABN device [9]. 
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Figure 2: (a) The high-tilt state, (b) The low-tilt state. 
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Figure 3: The post cross-section, L p x L p , and the cell cross-section, L c x L c , 
where L c = 2L p . 




Figure 4: The model geometry with L c — 1, L p — 0.5 and H = 3. When 
translated to physical units, these correspond to L c = 1/im, L p = 500nm 
and H = 3/im respectively. 
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Table 1: The four sets of vertical orientations. 

The last case necessarily creates a planar defect on this face and is, therefore, 
excluded (the constraint of interior continuity disallows certain choices of the 
edge orientations, for details see [H]). Of the remaining two cases, we choose 
the first one since it is more symmetric and is, therefore, expected to have 
lower energy. Once the horizontal orientations on the top face are fixed as in 
Figure the horizontal orientations on the corresponding bottom edges 
are taken to coincide with those on the top. This then completely determines 
the edge orientations on the eight horizontal post edges in our model, which 
are kept fixed for simplicity. 

We then consider the edge orientations on the four vertical post edges. 
We recall that tangent boundary conditions imply that n is either oriented 
upwards or downwards on these edges. Labelling the four edges by i = 
{1, 2, 3, 4} as in Figure (171) and the corresponding edge orientation by e^, there 
are four distinct cases up to symmetry and the sign of n. These cases are 
enumerated in Table (1). The first case, T, corresponds to n being oriented 
upwards on all four edges. For the second case, Pi, n is oriented upwards on 
three vertical edges and downwards on the remaining fourth edge whereas for 
the cases P2 and P3, n is oriented upwards on two edges and downwards on 
the remaining two. P2 and P3 are distinguished by the fact that changes 
sign on a pair of opposite faces in P2 whereas changes sign on all four 
vertical faces in P 3 . 

In addition to the edge orientations, the topology of n is also determined 
by its behaviour on the post faces and around the post vertices. We take 
an arbitrary path connecting a pair of adjacent post edges, on a post face. 
The tangent boundary conditions imply that as we move along this path, n 
rotates in the plane of the corresponding face. We only consider topologies 
whereby n undergoes minimal rotation between pairs of adjacent post edges 
(which is ±~ or a quarter-turn). In the terminology of [TT], we say that these 
topologies have zero kink numbers. 

The fixed horizontal edge orientations and the constraint of zero kink 
numbers allows us to qualitatively predict the behaviour of n on the bottom 
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Figure 5: (a) The most symmetric choice, (b) Admissible choice of lower 
symmetry. (c) Necessarily contains a planar defect. 
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Figure 6: (a) The profile on the bottom substrate, (b) The profile on the top 
face of the post. 



9 



substrate and the top of the post, as shown in Figure (jBk) and Figure (jHb)- 
On the bottom substrate, our choice of the horizontal orientations implies 
that n splits at one of the vertices (the bottom left vertex), follows the post 
edges and then rejoins at the top right vertex, aligning along the square 
cross-sectional diagonal. The net rotation between any two adjacent edges 
is just ±|, as is evident from Figure (jOk)- Similarly, on the top of the 
post (refer to Figure (JSb)), we have a fractional source defect at one of 
the vertices (the bottom left vertex), accompanied by a sink defect at the 
diagonally opposite (top right) vertex and any continuous n aligns along the 
corresponding diagonal. 

The sharp post vertices can be treated as fractional point defects in our 
model; the corresponding strength is measured in terms of the neighbour- 
ing distortion around the vertex. We only consider topologies where n has 
minimal distortion around post vertices or equivalently, these vertices are 
fractional defects of minimal degree (the solid angle subtended by n as it 
varies around a post vertex is just ±|). In the terminology of [IT] , these 
topologies have minimal trapped areas. 

This then defines four distinct topologies, {T, Pi, P 2 , P3}, all of which have 
fixed horizontal orientations as in Figure zero kink numbers and mini- 
mal trapped areas but are distinguished by their vertical orientations given 
in Table (1). The topologies are labelled by their vertical orientations for 
simplicity, since the vertical orientations are the only topological parameters 
in our model. 



4 Numerical Modelling. 

We now investigate the existence of equilibrium stable configurations with the 
four topologies cited in Section (3). To this end, we use the finite-element 
method to numerically solve the Euler-Lagrange equations PJ associated 
with the Oseen- Frank energy (T5]) [JO]. The numerical modelling is carried out 
in Femlab - a commercial partial differential equation solver based on the 
finite-element method [15J. 

(n, Vn) = A (^-w (n, Vn)j (4) 

where w (n, Vn) - the energy density has been defined in ([3]), nj is the i-th 
component of the unit-vector field and rtjj is the j-th partial derivative of 

We first construct trial configurations for each of the four topologies. 
Here, it is simpler to specify an unnormalized vector field N = (N x , N y , N z ). 
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The unit-vector field n is then given by 



N 

The fixed horizontal orientations allow us to have a single prescription 
for N x and N y in the four cases. 

N x = j mi, ( ™\\ x 1 ' " 



where L p is the post cross-sectional parameter and H is the cell height. 

For < z < /i, up to the height of the post, N z clearly has different forms 
to account for the different sets of vertical orientations. 

'zx(h-z), T, 
{z{h-z))x (l + cos (g) +«*(=;)), Pi, 
(z(h - z)) X cos fgV P 2 , 

(^(/i - z)) X COS ^gj X COS (g) P 3 . 



(7) 



For h < z < H (from the top of the post to the top substrate), N z is given 
by 

»-iFk ® 

for all topologies. (Note that for L c = 2L P and our domain - ^- < < 
|,0 < z < if, the only points where N vanishes are given by the post 
vertices, which are necessarily singularities because of the tangent boundary 
conditions.) 

Given equations (jSJ), © and (jSJ), one can easily check that the represen- 
tative (j3J) satisfies the boundary conditions with the correct topology and is 
continuous away from the sharp post vertices. We use these trial configura- 
tions as initial conditions for the numerical solver. We numerically solve the 
Euler-Lagrange equations on a variable mesh that allows for greater resolu- 
tion near the sharp post features. The numerical solver respects the bound- 
ary conditions and topology so that the numerical solutions have the same 
topology as the initial condition. 

We show the solution profiles with h = L c in Figures (jHJ), (jSJ) and ( fill) . We 
plot the profile around the rectangular post for {T, Pi, P 2 , P3} in Figures 
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- d). It is evident that the profiles are continuous and regular around the 
post with no point or line singularities. In particular, there are no defects 
along the leading and trailing edges, as suggested by previous modelling 
[9]. In Figures (Eh - d), we describe the four profiles on the post faces in 
greater detail. Figure ([9k ) clearly shows that the T profile tilts upwards at 
all points as we move across the post faces and aligns along the post diagonal. 
Figures (|Hb-d) correspond to the Pi, P 2 and P 3 profiles respectively. Here, 
we note that whenever the vertical orientation changes sign between a pair of 
successive vertical edges on a post face, there is necessarily an intermediate 
planar region. This planar region is typically a curve across the face, between 
the two edges in question. We can clearly see these planar regions on a pair 
of adjacent faces in the Pi profile, on a pair of opposite faces in the P 2 profile 
and on all four vertical faces in the P3 profile. 

In Figure ( ITTi) . we plot the solutions on a periodically extended domain, 
comprising of a pair of posts, along the cross-section y = L c /4 as shown 
in Figure (TlO|) . The periodic boundary conditions imply that the profiles 
are simply repeated between neighbouring posts. In particular, the solution 
topology is the same for every post. This is illustrated in Figure (FiTh ). where 
the T profile tilts upwards between a pair of neighbouring posts, whilst inter- 
polating between the fixed orientations on the successive vertical edges. Sim- 
ilarly, the planar regions observed in Figures ([9b-d) extend between neigh- 
bouring posts in Figures ffTTb-d). by virtue of the periodic boundary condi- 
tions. Above the post, for h < z < H, all profiles tilt upwards and finally 
become homeotropic on the top substrate. These profiles, though discussed 
for a fixed post height h = L c , are characteristic of the solution profiles for 
all h. 

Next we look at how the solution energy varies with scaled post height, 
h/L c , for two different sets of elastic constants in Figures (|T2l) and (|TS1) . 
The second set, with equal elastic constants, corresponds to the widely used 
one-constant approximation [U [TO]. It is evident from Figures (fl2l) and 
(fl3l) that the qualitative trends for the two sets of elastic constants are the 
same. However, the one-constant approximation in Figure ffl3|) is found to be 
computationally less demanding than the unequal constant case and allows 
for greater numerical resolution. 

We first note that the T solutions have the smallest energy whereas the P3 
solutions have the highest (there is a slight crossover between the P 2 and P3 
solutions in Figure (fl2|) but the free energy difference is negligible). This is 
consistent with experimental observations which show that the liquid crystal 
always relaxes into the high-tilt state, when cooled down from the isotropic 
state [S]. 

These energy trends can be anticipated and explained on topological 
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Figure 8: (a) The T profile, (b) The P x profile, (c) The P 2 profile, (d) The 
P 3 profile. 
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Figure 9: The solution profiles on the faces, (a) T (b) P\ (c) P2 (d) P3 
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grounds. As is evident from Figures (jHfci) and (fTTa ). the T solution has mini- 
mal distortion consistent with the boundary conditions. The P3 solution, on 
the other hand, has the most strongly distorted profile because of the planar 
regions on every vertical post face and between every pair of neighbouring 
posts. These planar regions arise from the continuous interpolation between 
opposite vertical orientations; these vertical orientations, in turn, are part of 
the topological characterization of these solution profiles. The minimal and 
maximal distortion, in the T and P3 solutions respectively, then qualitatively 
explain the energy trends. 

We look at the free energy differences in Figure (1141) . Here, we plot two 
quantities 

e i = (9) 



E 



T 



and 



HiT 

for the elastic constants in Figure (1T21 . E Pl is typically the smallest of 
{E Pl , Ep 2 , Ep 3 } and E P , A is the largest. It is obvious from ffT4"l) that both 
€\ and 63 can be qualitatively described as piecewise linear functions of 
with the gradient changing as a function of post height. The gradient is 
smallest in the interval -p- £ [0.5,0.9], where both e\ and 63 are small and 
grow very slowly. We refer to this interval, h £ [0.5 x L c ,0.9 x L c ], as the 
plateau region in the remaining discussion. For post heights h > L c , the 
gradient of both curves increases appreciably so that ei and €3 grow relatively 
quickly with post height. This change in behaviour as a function of post 
height can be explained by Figure ( |T2l) . Et, the energy of the tilted solution, 
typically reaches a shallow local maximum and then saturates as the post 
height increases. On the other hand, both Ep x and Ep 3 grow linearly with 
post height. These differences in the energy trends manifest in e\ and e3, 
given by (J9j) and ( JTU1) . leading to the features cited above. 

To understand what these numerical calculations mean in the context 
of the actual device, we first note that the tilted and planar states are 
observed experimentally for post heights in the range 0.6 microns to 1.2 
microns [private communication]. This bistable region corresponds to h £ 
[0.6 x L c , 1.2 x L c ] in our model. Further, as the posts become taller, the 
lower energy tilted state is more easily observed than the planar state and 
it becomes increasing difficult to switch from the tilted state to the planar 
state p. Our modelling shows a plateau region for £ [0.5,0.9], where 
both €\ and 63 are small and do not change significantly. This plateau region 
suggests that we might observe the topologically distinct T and Pi states, 
i = 1,2,3, for £ [0.5,0.9]. Our numerical range and the experimentally 
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found bistable region are in qualitative agreement. Secondly, the fact that 
ei and e 3 grow quickly as a function of y 5 - indicates that the Pj states be- 
come energetically far more expensive than the T states as the post height 
increases. This, in turn, qualitatively explains the experimental observations 
and switching characteristics for taller posts. A more detailed analysis of the 
liquid crystal energy would include surface effects and scalar order parame- 
ters mini]. 

5 Conclusion. 

The methods and results presented in this paper illustrate how topology can 
be exploited in order to find different static liquid crystal configurations in 
prototype device geometries. For the specific case of the PABN geometry, we 
have identified four simple, low-energy topologies. This topological informa- 
tion is then used to construct suitable initial conditions for a finite-element 
numerical algorithm, yielding four distinct classes of numerical solutions. It 
would be difficult to numerically find these different solutions without the 
topological insight. Our modelling suggests that the tilted and planar states 
in the PABN geometry are topologically distinct, leading to a topological 
mechanism for the observed bistability. (For example, the T topology de- 
fined in Section (3) is a good candidate for the topology of the tilted state 
whereas the P, topologies serve as good candidates for the planar state.) Sec- 
ondly, these numerical solutions have comparable free energies over a certain 
range of post heights, indicating a limited bistability region that is again 
commensurate with experimental observations. 

The methods outlined in this paper can be extended to more general 
polyhedral geometries. We can use topology to organize the space of ad- 
missible configurations, understand the structure of the stable configurations 
and their energetics. This will be investigated further in subsequent work. 
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Figure 10: A periodically extended domain containing two posts. This is 
the bottom cross-section, with the filled regions corresponding to the post 
bases. The bold line indicates the cross-section y = L c / 4, which is plotted in 
Figure (HH). 
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Figure 11: (a) The T profile, (b) The P\ profile, (c) The P2 profile, (d) 
The P3 profile. These plots are the projection of the solution profiles onto 
the plane y = L c /4 i.e. only the x and z components of the profiles are 
shown here. The dots, in particular, correspond to points where the solution 
is oriented in the ^-direction. Note that we only plot up to height z = 3h/2 
because the profiles have similar structure in the remainder of the cell ^ < 
z < H, where h is the post height and H is the cell height in Figure (jlj). 
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